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Abstract 

In this paper we study path integral for a single spinless particle on a star graph with 
N edges, whose vertex is known to be described by U{N) family of boundary conditions. 
After carefully studying the free particle case, both at the critical and off-critical levels, 
we propose a new path integral formulation that correctly captures all the scale-invariant 
subfamily of boundary conditions realized at fixed points of boundary renormalization 
group flow. Our proposal is based on the folding trick, which maps a scalar-valued 
wave function on star graph to an iV-component vector-valued wave function on half- 
line. All the parameters of scale-invariant subfamily of boundary conditions are encoded 
into the momentum independent weight factors, which appear to be associated with the 
two distinct path classes on half- line that form the cyclic group Z2. We show that, when 
bulk interactions are edge-independent, these weight factors are generally given by an 
A^-dimensional unitary representation of Z2. Generalization to momentum dependent 
weight factors and applications to worldline formalism are briefly discussed. 
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1 Introduction 



Over the last one and a half decade, quantum graphs have been attracted much atten- 
tion as the simplest low-energy effective theory for mesoscopic networks of one-dimensional 
quantum wires, whose configuration space is effectively modeled by a graph (see for a brief 
survey [1] and references therein). In general, a quantum graph is a non-differentiable man- 
ifold and has singularities at junction points (vertices) of graph edges. A standard lore for 
taming such singularity is first to require Kirchhoff 's law of conserved currents (Noether 
currents) at the vertex and then to reduce to the problem of boundary conditions for wave 
functions/quantum fields defined in the bulk. 

Although Kirchhoff 's law of conserved currents has been vastly studied in the literature 
in the operator formalism, its path integral description is not yet fully understood. The 
problem is that boundary conditions provided by Kirchhoff 's law of conserved currents will 
be characterized by a large number of parameters, however, there is no a priori prescription 
to include these parameters into the path integral weight exp(i5) nor into the path integral 
measure Dx. Although there exist several attempts, most of the literature concern path 
integral description for a system defined on a half-line with the Robin boundary condi- 
tion [SHE], which is characterized by only one real parameter, or a system on a line with 
a Dirac (5-potential [91-115) . which is characterized by a single real coupling constant, or a 
system on a circle with a single scale-independent point interaction [71ll61[T7|. which is char- 
acterized by two real parameters, or a system on a circle with a single U{2) family of point 
interactions [181I19| . which is characterized by four real parameters. Path integral approach 
to more complicated quantum graph that contains more boundary condition parameters is 
yet mysterious, even for the case where there is no interaction in the bulk. In this paper 
we would like to tackle with this problem. To get the feel of our problem, in the first part 
of this paper we will concentrate ourselves to non-relativistic quantum mechanics for a free 
spinless particle on a star graph with edges, whose vertex is known to be described by 
U{N) family of boundary conditions and hence characterized by A^^ real parameters. Since 
the Schrodinger equation on star graph is a well-defined problem, the most unambiguous 
way to find the path integral formulation on star graph will be as follows: (i) first solve 
the Schrodinger equation with U{N) family of boundary conditions, (ii) and then find the 
complete orthonormal set of energy eigenfunctions, (iii) and then evaluate the Feynman 
kernel in terms of the complete orthonormal set, (iv) and finally rewrite the kernel into the 
path integral representation. With this approach we will show that all the A^'^ parameters 
are naturally encoded into the weight factors of scattering matrix (S-matrix) associated 
with the two distinct path classes that form the cyclic group of order 2 (Z2). In the second 
part of this paper we will proceed to include edge-independent bulk interactions and show 
that, when the theory lies on a fixed point of boundary renormalization group (RG) flow, 
the weight factors are generally given by an A^-dimensional unitary representation of Z2. 

The rest of this paper is organized as follows. In Section [2] we first recall the complete 
orthonormal set of energy eigenfunctions of the free Hamiltonian on star graph and then 
compute the Feynman kernel. By taking the continuum limit of time-slicing representation 
of the kernel, we find that the path integral description of U (N) family of boundary con- 
ditions is naturally formulated into the phase space path integral rather than configuration 
space path integral with exponential weighting factor exp{iSfree[x{t),p{t)]), where Sf^ee is 
the action for a free particle in the Hamiltonian formulation, and another weighting factor 
of S-matrix §(p) that depends on the initial momentum p = p{0). The resultant path in- 
tegral becomes a summation over two distinct path classes categorized by the number of 
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(a) Unfolding picture. (b) Folding picture. 



Figure 1: Star graph with N edges. In the folding picture wave function of the system is 
given by an A^-component vector of L^(M+)-functions. 

interactions (n = or 1) for which a classical particle experiences at the vertex. Contribu- 
tions from the path class of n = 1 turn out to be weighted by the S-matrix. In Section [3] 
we discuss the critical case; that is, the case of boundary conditions realized at the fixed 
points of boundary RG flow. In this critical case the S-matrix becomes constant matrix 
such that the phase space path integral can be simply reduced to the ordinary configura- 
tion space path integral. In Section U] we generalize to an interacting case by adding an 
edge-independent bulk interaction. We show that, at the fixed point of boundary RG flow, 
the weight factors are generally given by an A^-dimensional unitary representation of the 
cyclic group Z2, which is the main result of this paper. Section [S] is devoted to conclu- 
sions, discussions for remaining issues and a simple application to worldline formalism on 
star graph. Computational details for the proof of orthonormality and completeness of the 
eigenfunctions of the free Hamiltonian on star graph are presented in Appendix \^ 



2 Free particle 

To begin with, let us first recall non-relativistic quantum mechanics for a free particle on a 
star graph with edges. The star graph T is the one-point union of N positive half-lines 
M.f = {xi G M I Xj > 0} (i = 1, • • • , A^) such that the Hilbert space T-L^T) is just given by 
the direct sum V. = 0,^ii^(M+). Since L^{Rf) is equivalent to L'^{Rp for any the 
Hilbert space can also be written as follows: 

^ = l2(M+)®C^. (2.1) 

It should be emphasized that the concept behind the change of notation from L^(M^) 
to L2(k+) cN so-called folding trick [20H23], which maps a single "scalar-valued" 

wave function on star graph to an A^-component vector-valued wave function on half-line; 
see Figure [H In the folding picture the wave function (element oi T-L) is then given by 

tf{x) = {Mx),--- ,i^N{x)f, (2.2) 

where tpj{x) (j = 1, • • • , N) corresponds to the wave function on the jth edge with x being 
the distance from the vertex. T stands for transposition. The time-independent Schrodinger 
equation for a free particle now becomes the following matrix equation 

^iree^{x) = Elfix), X / 0, (2.3) 
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where the bulk Hamiltonian is given by the N x N diagonal matrix 

Hfree := diag(i?free, • • • , -f^free) = -^free ^ I, -f^free : = 



(2.4) 



with I being the N x N unit matrix. In this paper we will work in the units where h = 
2m = 1. As is well-known, Kirchhoff's law of probability currents (or self-adjoint extension 
of Hfree) leads to the following U{N) family of boundary conditions at the vertex: 



(I - U)^(O) - iLo(I + U)^'(O) = 0, U G U{N), 



(2.5) 



where prime (') indicates the derivative with respect to x. Lq is an arbitrary length pa- 
rameter (or renormalization scale j24) ) which is inevitably introduced to adjust the length 
dimension of the equation ()2.5p . In this paper we will use black board symbols (A, B, C, 
etc.) to denote matrices. For the following discussions it is convenient to parameterize the 
matrix U € U{N) into the following spectral decomposition form 



U 



N 
a=l 



(2.6) 



where e*"" (0 < a" < 27r) is the ath eigenvalue of U and ^'^ = ■ ■ ■ is the corre- 

sponding ath eigenvector that fulfills the eigenvalue equation U^'^ = e*" the orthonor- 
mality (I*'*)"'" ■ ^ = 5^'^ and the completeness '^a=i^'^ ' i^'^)^ — ^- With these eigenvectors 
P° becomes the hermitian projection operator satisfying "^^^i P"^ = I, P'^P^ = 5"*P* and 
(P")t = p'^j where dagger (^) stands for the hermitian conjugate. By substituting ()2.6p into 
the equation ()2.5p . the boundary condition is diagonalized and boils down to the following 
independent equations 



(Ot. v(o) + LV(o) 



0, := Ln cot 



a 



(2.7) 



Now it is easy to find the solution to the Schrodinger equation (|2.3p with the boundary 
condition (|2.7p . Let us first consider positive energy (E > 0) solutions. For given momentum 
p = \fE > 0, the continuum spectrum exhibits A^-fold degeneracy and there exist distinct 
eigenfunctions which are orthogonal to each other yet have the same energy eigenvalue 
E = p'^. One of the most convenient choices of such solutions is as follows: 



'il)"'{x;p) = [lexp(— ipx) -|- S(p) exp(-l-ipx)],^" 



exp(— ipx) -|- 



ipL" - 1 



exp(-l-zpx) 



ipL" + 1 

where S(p) is the one-particle S-matrix given by 

N 



C 0<p<oo, 



^N, (2. 



Hp) = E 



a=l 



ipL"- + 1 



exp 



N 



', 2arccot(pL° 



a=l 



(2.9) 



whose iz-component gives the reflection coefficient for a particle residing on the ith edge and 
ij-component (i ^ j) the transmission coefficient for a particle traveling from the jth edge 
to the ith edge. Note that the S-matrix satisfies the unitarity, E>^{p)S{p) = E>{p)§'^{p) = I, 
and the hermitian analyticity, S^(p) = S{—p). Consequently, ^°(x;p) and il^°'{x;—p) are 
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Figure 2: Integration contour C and simple poles of one-particle S-matrix S(p) on the 
complex p-plain. 



just related by the unitary transformation 'ip^{x;p) = §(p)^"(x; — p). Note also that any 
solution to the Schrodinger equation ()2.3p with the energy eigenvalue E = is given by 
the linear combination oi 'il)"'{x]p), a = 1, - ■ ■ ,N. 

Let us next consider negative energy {E < 0) solutions, or bound states, which are 
characterized by the S-matrix pole lying on the positive imaginary p-axis. As discussed 
in |24) . there is a one-to-one correspondence between the eigenvalues of U and the negative 
energy spectrum: If the eigenphase a° lies on the range < a" < vr (or < L"^ < oo), 
there exists a normalizable bound state with the energy E = — 1/(L'^)^. If the eigenphase 
a" lies on the range tt < a"" < 2it (or — oo < < 0), on the other hand, there exists a 
non-normalizable antibound state with the energy E = — 1/(L'')^. The normalizable bound 
state solution that corresponds to the eigenvalue e*"" is found to be of the form 




exp 



X 

1? 



C, < L° < oo. 



(2.10) 



The set of these eigenfunctions (|2.8|) and (j2.1Up provides the complete orthonormal basis of 
the system; that is, they satisfy the orthonormality 



dx [ip^ixip)]^ ■^ix;q) = 2TT5''^6{p-q), 



dx [rsix]]^ ■i?s{x) = 5'^\ 



dx[nx;p)y ■^Bix)= / dx[rBix)V ■^\x;p) = 0, 
Jo 



(2.11a) 
(2.11b) 
(2.11c) 



and the completeness 

N 



J2 r^r{x;p)-[r{y;p)y + Yl rB{x)-[rB{y)]^ = S{x-y)l. (2.12) 

a=l •'^ 0<L''<oo 

(For the computational details, see Appendix|Xl) By using this complete orthonormal basis 
the Feynman kernel 'K[x,y;T) = (x|e~*'^f'^'=°^|y) can be expressed as the following N x N 
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matrix 



a=l ^ 



0<L"<oo 



exp 



(2.13) 



whose ij-component gives the transition amphtude for a particle traveling from {y,j) to 
(x, i) in time interval T, where the notation (x, i) stands for the position located at the 
distance x from the vertex on the ith edge. As easily checked the Feynman kernel ()2.13p 
satisfies the following equations: 



(i) Schrodinger equation [idx^ — Mfj.cc{—idx)]^{x , y; T) = 0; 



(ii) Initial condition 

(iii) Composition rule 

(iv) Unitarity 

(v) Boundary condition 



K(x,y;0) =<5(x-y)I; 

poo 

/ dz ]K(x, z- Ti)K{z, y- T^) = K{x, y; Ti + Ta); 
Jo 

K\x,y;T)=Kiy,x;-T); 

(I - U)K{0, y; T) - iLo(I + U)(a,K)(0, y; T) = 



0. 



(2.14a) 
(2.14b) 
(2.14c) 
(2.14d) 
(2.14e) 



By substituting the solutions (j2.8p (j2.10p into (j2.13p the Feynman kernel can be written as 
the following momentum integral 



K(x,y;r) 



dp 



,27r 



dp 



c 



2tt 



d> -^\K{x,y;T\p), 



(2.15) 



where IK is the momentum representation of the kernel and given by 

'x - y" 



K{x, y;T \ p) = I exp iT 



P 



T 



P 



+ S(p) exp ( iT 



P 



x + y 



P 



(2.16) 



The integration contour C is chosen to enclose all bound state poles on the complex p- 
plain in clockwise direction; see Figure [2j The contour integral picks up the bound state 
contributions to the Feynman kernel such that it will vanish if the S-matrix has no pole on 
the positive imaginary p-axis. 

The physical interpretation of (j2.16p is now clear: the exponent T[p{{x =p y)/T) — p"^], 
which is invariant under the parity {x,y,p) i— )• {—x,—y,—p), is nothing but the classical 
action for a free particle propagating along the classical trajectory Xci{t) = ±y + {{x^y)/T)t 
with fixed conserved momentum Pci{t) = p, where for Xci{t) = —y + ((x + y)/T)t < it 
should be regarded as a mirror image Xc\{t) = y — {{x + y)/T)t with flipped momentum 
Pc\{t) = —P'-, see Figure El and m 

By using the kernel ()2.15p time-evolution of any state "0 S L^(M+) (g) that satisfies 
the boundary condition (j2.5p is given by 



V^(x,r) 



dy 



dp 



+ 



c 



dp 
2^ 



K(x,y;r I p)V'(y). 



(2.17) 
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Figure 3: Classical worldlines for a free particle on a star graph in the unfolding picture. 



Now we wish to represent ]K(x, y;T \ p) as a functional integral over all possible paths 
fluctuated around the classical trajectories. To this end, let to = < ti < • • • < < 
tn+i = T (ti = £5t with 5t := T/{n + 1)) be the equidistant time sequence. With this time 
sequence the exponential weighting factor exp(iT[p((x =F y)/T) — p^]) can be expressed as 
the following time-slicing representation: 



exp 

n 

n 



p 



x=fy 



P 



.1=1 



dxf 



exp I 



^6t 



Pe 



V St 



p} 



Xn + l=X 



, (2.18) 



xo=±y; po=P 



which is of course invariant under the parity {x£+i, Xi,pi) i— )• {—xe+i,—Xi,—pi). It should 
be emphasized that the integration range of x-coordinate is enlarged to the whole line M 
such that the path integration must be performed over paths including those depicted in 
Figure 4(c) and |4(d)] Now, by taking the limit n — t- oo (keeping T fixed) we get the following 
phase space path integral over the path {x,p) : [0, T] — )• M x M with the boundary conditions 
x(0) = ±y, x(T) = X and p{0) = p: 



Kix,y;T\p)=I 



x{T)=x 

Vx{t) 

x{0)=y Jp{0)=P 
f-x{T)=x 
+ §{p) / Vx{t) 
Jx{0)=-y 



Vp{t) exp(zS'frec[a:(t),p(t)]) 



p(0)=p 



Vp{t) exp(i5frcc[2;(t),p(t)]). 



(2.19) 



where S{T:e(.[x{t),p{t)] is the action for a free particle in the Hamiltonian formulation 



Slree[x{t),p{t)] 



dt [p{t)x{t) - H{,ee{p{t))] , //frce(p(i)) = p\t), 



(2.20) 



with dot ( ' ) indicating the derivative with respect to t. (Recall that we are working in the 
units h = 2m = 1). The integration measures are given by 



r roo 

Vx{t) := lim / dx. 

/ Vp{t) := lim TT / 

/ n—^oo / 



2tt 



(2.21a) 
(2.21b) 
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(a) Direct path (weight factor: Sij). (b) Reflected path (weight factor: Sij). 




(c) Typical quantum fluctuation of (a) (d) Typical quantum fluctuation of (b) 
(weight factor: Sij). (weight factor: Sij). 



Figure 4: Classical worldlines and its typical quantum fluctuations for a free particle on a 
star graph in the folding picture. Time is flowing along the vertical direction. 



Notice that the total numbers of dxg and dp£ in the measures are the same as opposed to 
the standard phase space path integral, where Dp contains one more dpi compared to Dx. 

We emphasize that, although rewriting the kernel ()2.13p into the path integral ()2.19p 
is an almost trivial exercise, equation (j2.19p is very suggestive in a sense that all the N'^ 
parameters which characterize the vertex are brought by the unitary weight factors {I, 
associated with the two distinct path classes on the half-line. This fact strongly encourages 
us to understand the unitary matrices {I, S(p)} as a unitary representation of a certain 
discrete group that all possible paths on M"^ form, just like the standard weight factors in a 
non-simply connected space |25H32j . This naive expectation will be more transparent and 
turn out to be true when the S-matrix becomes independent of the momentum p, which is 
the main subject of this paper and will be discussed in detail below. 



3 Free particle at criticality 

A remarkable simplicity occurs when the S-matrix does not depend on p, which can be 
realized as fixed points of boundary RG flow. Such boundary RG flow has first been 
studied in ^33„34J (but with different interpretations of UV and IR fixed points) and later 
in [21] . Search for universality classes of boundary conditions is very important not only for 
theoretical interests but also for physical applications to transport phenomena of Tomonaga- 
Luttinger (TL) liquids on quantum wire junctions. Enormous studies have been conducted 
to elucidate the rich fixed point structure of TL liquid junctions (see for example [351-I3H]). 
however, there is no universally accepted phase diagram in the community thus far. A 
geometric approach to scale-invariant boundary conditions which we will present in the 
next Section may shed new light on this subject. We note that scale-invariant subfamily of 
boundary conditions have also been discussed in the literature |16 p i7| ,[39VI¥I] independently 
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of the renormalization group study. 

Now, let us move on to the analysis of path integral at criticality. As shown in 
boundary conditions at the fixed points are characterized by a hermitian unitary matrix 
satisfying the condition = I and parameterized by A^^ — N = N{N — 1) real parameters. 
There are 2^ distinct fixed points in the theory space, whose relevancy or irrelevancy are 
classified by the number of +l-eigenvalues (or — 1-eigenvalues). At those fixed points, the 
S-matrix becomes identical to U and the momentum integrations ()2.18p and ()2.15p are easily 
performed. The resultant kernel can be written as the following configuration space path 
integral over the path x : [0,r] — )• M with the boundary conditions x(0) = ±y, x{T) = x: 



K{x,y;T)=I 



+ U 



x{T)=x 

Vxit) exp[iS{rcc[xit)]) 

x{0)=y 
x{T)=x 

Vx{t) exp(iS'frcc[a;(t)]), 

x{0)=-y 



(3.1) 



where S{ree[x{t)] is the action for a free particle in the Lagrangian formulation and given by 

SircMt)] = [ dth\t). (3.2) 



The integration measure now becomes that for standard configuration space path integral 
(so-called Feynman measure) 



/ 



Vx(t) :-- 



lim ^ 



n 



dxf 



5t 



n + l 



(3.3) 



Notice that the boundary condition for the kernel ()3.ip becomes two independent equations 
(I-U)]K(0, y;T) =0 and (I+U)((9^1C)(0, y; T) = 0, which follow from the boundary condition 
and the orthonormality relations [(I±U)/2][(I=fU)/2] = and [{l±\J)/2f = (I±U)/2 
for a hermitian unitary matrix U; see also the discussion below (|4.8ap and ()4.8bp . 

Now it is easy to perform the path integration (|3.ip . or, equivalently, momentum inte- 
gration ()2.15p with ()2.16p . A straightforward computation gives 



]K(x,y;r) 



lexp 



y 



T 



-MJexp 



,t1 

4 



X + y 



(3.4) 



Notice that the exponent (T/4)((x =F y)/T)'^, which is again invariant under the par- 
ity {x,y) I— {—x,—y), is nothing but the classical action S'frecla^cK*)] for a free parti- 
cle traveling along the classical trajectory Xci(t) = ±y -|- {{x =F y)/T)t. The overall fac- 
tor l/\/47riT, on the other hand, is the standard contribution from the path integral 
/g(o)=g(T)=o ■^'^(^) ^"'^P(^'^free[^(*)])' which follows from the redefinition of the path x{t) = 
Xci{t) + lit), such that it can be regarded as "one-loop correction" [32]. As we will show 
in the next Section, N x N matrices {I, U} can be derived as an A^-dimensional unitary 
representation of the cyclic group Z2, whose origin is the geometry of half- line. In this sense 
quantum mechanics for a free particle on a star graph at a fixed point of boundary RG flow 
can be regarded as "classical" -|- "one-loop correction" -|- "Z2 factor". 

We note that when A^ = 1 we have only two choices for the hermitian unitary weight 
factor, U = -|-lorI[J = — 1, the former corresponds to the Neumann boundary (UV flxed 
points of boundary RG flow [23]) and the latter the Dirichlet boundary (IR fixed points of 
boundary RG flow [23]), which coincide with the well-known results. 
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(a) [even]: Class of loops that cross (b) [odd]; Class of loops that cross 
the origin even number of times. the origin odd number of times. 




(c) [odd] * [odd] — [even]: The first arrow follows from Z2 symmetry of the upper loop. 



Figure 5: Two distinct loop classes on M/(x ~ —x), where the base point x is chosen to 
exclude the origin x = 0. (The case x = falls into the problem of boundary conditions 
()4.8ap and (|4.8bp .) Composition of these loop classes satisfies the Z2 multiplication rules: 
[even] * [even] = [odd] * [odd] = [even], [even] * [odd] = [odd] * [even] = [odd], where * 
indicates composition of paths. 

4 Adding bulk interaction: A geometric perspective 

So far we have studied path integral formulation for a free particle on a star graph by 
using the time-slicing representation of the kernel obtained from the operator formalism. 
We have shown that all of the information about the boundary conditions can be encoded 
into the weight factors. From the path integral standpoint, however, it is more desirable 
to formulate the weight factors without using the eigenfunctions of Hamiltonian nor the 
boundary conditions for wave functions. In this Section we will show that, without referring 
to the boundary conditions (j2.5p . the momentum independent weight factors {I,U} can be 
derived as an A^-dimensional unitary representation of the cyclic group Z2 = {/, P}, where 
/ (identity) and P (parity) satisfy the following multiplication rules 

lo I = Po P = I, IoP = PoI = P. (4.1) 

Here circle (o) stands for abstract group multiplication. / and P act on the coordinate 
space as usual, I : x ^ x and P : x ^ —x. We will see that / corresponds to the unit 
matrix I and P the hermitian unitary matrix U. 

To this end, it is worthwhile to remember first the standard argument of weight factors 
in path integral, which can be summarized as the Schulman-Laidlaw-DeWitt theorem |25l 
[25] (see also |27fl32] ). In a non-simply connected space Ai, the total amplitude for a 
given transition is obtained by a linear combination of partial amplitudes, where each 
partial amplitude is constructed from configuration space path integral on the universal 
covering space of M over paths that belong to the same homotopy class. The Schulman- 
Laidlaw-DeWitt theorem states that the coefficients in the linear combination must be a 
one-dimensional unitary representation of the fundamental group vri(A^)Il] Although the 

^Horvathy et al. discussed a generalized case where, if wave function becomes vector- valued due to some 
internal symmetries, weight factor could become a higher-dimensional unitary representation of 7ri(A^) [30| . 
This is very similar to our situation although there is no genuine internal symmetry in our problem. 
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half-line is a simply-connected space and its fundamental group is trivial, as discussed in the 
previous Sections and is widely known under the name of method of images, path integral 
on M"^ should be lifted to that on the whole line M with an identification x ~ — x and given 
by a linear combination of partial amplitudes associated with the two distinct path classes 
on M/(x ~ —x), one is a class of paths that cross the origin even number of times and 
the other odd number of times [3]. It is intuitively clear that, for fixed initial and final 
positions {x{T),x{0)), these two path classes form the cyclic group of order 2 under the 
composition of paths, and are distinguished by Z2-action on the initial position x(0) = y 
{y > 0); see Figure [5] for the case of loops (x(T),x(0)) = (x,±x). Now that space of all 
possible paths we should integrate over is equipped with the group structure Z2, along the 
same philosophy of the Schulman-Laidlaw-DeWitt theorem, the weight factors in the linear 
combinations must fall into the one-dimensional unitary representation of Z2. There are 
only two possibilities of such unitary representations, {+1,+!} (trivial representation) or 
{+1,-1} (defining representation), the former gives the weight factors for the Neumann 
boundary and the latter for the Dirichlet boundary. Similar arguments are given in [43J by 
regarding the half-line as an orbifold M/Z2 from the beginning. 

Now we wish to extend the Schulman-Laidlaw-DeWitt theorem to the path integral on 
a star graph. An important point to note is that, in the folding picture, the configuration 
space is but the Hilbert space becomes Ti = L^(M+)(8)C^ such that the Feynman kernel 
IC that acts on Ti must be an x matrix function. Furthermore, if the bulk interaction 
is edge-independent and Hamiltonian M is given hy M = H <^ I, which is the case of free 
particle, H and boundary conditions must be simultaneously diagonalizable such that the 
kernel should be factorized as IK = (A^ x A^ constant matrix) x (scalar kernel). Thus we can 
start with the following ansatz 

K{x,y;T)=W{I)K^{x,y;T) + W{P)K^{x,-y;T), x,y > 0, (4.2) 

where W(/) and W(P) are N x N constant matrices to be determined later. K^{x,y;T) = 
{x\e~^^'^\y) is the scalar kernel defined on the whole line M. and assumed to satisfy the 
following conditions: 

(i) Schrodinger equation [idx — H{x, —idx)]KR{x,y;T) = 0; (4.3a) 

(ii) Initial condition K^{x,y;0) = 5{x — y); (4.3b) 

/OD 

dz K]^{x,z;Ti)Kt>{z,y;T2) = Kf>{x,y;Ti + T2); (4.3c) 
-00 

(iv) Unitarity K^{x,y;T) = K^iy, x; -T); (4.3d) 

(v) Z2 symmetry Ku{x,y]T) = K^i-x, -y;T), (4.3e) 

where asterisk (*) indicates complex conjugate. Equation ()4.3ap assures the total amplitude 
()4.2p to satisfy the Schrodinger equation [idx^ — M{x, —idx)]^{x,y;T) = 0. We note that 
Z2 symmetry (|4.3e|) will be justified for parity invariant bulk Hamiltonian H. Construction 
of scalar kernel satisfying the conditions (|4.3aP " (|4.3ep is a standard problem so that we 
assume is given by an ordinary configuration space path integral. 

With these assumptions, the total amplitude IK will satisfy the initial condition (|2.14bp . 
composition rule ()2.14cp and unitarity ()2.14dp if and only if the weight factors form an 
A^-dimensional unitary representation of the cyclic group Z2. Let us prove this statement 
by imposing these conditions to the kernel IK. 
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• Initial condition. Let us first impose the initial condition ()2.14bp to K. A straightfor- 
ward calculation gives lC(x, y; 0) = W{I)5{x — y) + W{P)5{x + y), where we have used 
the assumption ()4.3bp . Since 5{x + y) = for x,y > 0, the initial condition ()2.14bp 
leads to the following constraint 

W(I) = I. (4.4) 

• Composition rule. Let us next impose the composition rule ()2.14cp to K. By substi- 
tuting the ansatz the left hand side of ()2.14cp becomes /g°°dzK(x, z; Ti)]K(z, y; T2) = 
(W(/)W(I) /~dz+W(P)W(P) f_^dz)K^{x, z- Ti)Kuiz, y; r2) + (W(/)W(P) J^dz+ 

W(P)W(/) f^^dz)KK{x,z;Ti)Kf>{z,-y;T2), where we have used g^H)- In order to 
realize the right hand side of (|2.14cp the weight factors must obey the following mul- 
tiplication rules 

W(/)W(/) = W(P)W(P) = W(/), (4.5a) 
W(/)W(P) = W(P)W(I) = W(P), (4.5b) 

which are nothing but those for Z2; see ()4.ip . Thus the weight factors {W(/), W(P)} 
must be an X matrix representation of the cyclic group Z2. 

• Unitarity. Let us finally impose the unitarity condition ()2.14dp to K. By taking 
the hermitian conjugate of the ansatz we get K\x,y;T) = W\l)Ku{y,x; -T) + 
W^P)KK{y, —x; —T), where we have used ()4.3dp and (|4.3ep . Thus the unitarity 
leads to the following constraints 

W\l) = W(/), W\P) = W(P). (4.6) 

Hence the weight factors {W(/),W(P)} must be hermitian matrices. It follows im- 
mediately from ()4.4p . ()4.5ap and ()4.6p that the weight factors must also be unitary, 
W-^(/) = Wt(/) and W^iP) = W^P). 

Collecting the above results we conclude that W must be an A^-dimensional unitary repre- 
sentation of the cyclic group Z2; that is, it must be a map W : Z2 — )• U{N). Such unitary 
representation is given by 

W(/) = I, W(P) = U, (4.7) 

where U G U{N) is a hermitian unitary matrix that satisfies = I. Different choice of U 
corresponds to different quantization. 

We emphasize that the weight factors are determined without referring to the boundary 
conditions. This is very important because the path integral presented in the previous Sec- 
tion is based on the boundary conditions provided by the self-adjoint extension of free Hamil- 
tonian, while the results (|4.7p do not depend on any specific Hamiltonian (except for the par- 
ity invariance). Indeed, boundary conditions for IK can be derived as a consequence of ()4.7p : 
If the scalar kernel is continuous and smooth at the origin, ^15(0+, y; T) = A'k(0_, y; T) and 
(5a;-R'K)(0+, y; T) = {dxK^){Q-,y:,T), which will be justified when there is no singularity 
at X = in the parity invariant bulk Hamiltonian, we immediately obtain the following 
boundary conditions for the total Feynman kernel: 

(I - U)IK(0+, y; T) = (I - U)(I + V)K^{Q+,y- T) = 0, (4.8a) 
(I + U) {dJL) (0+ , y ; T) = (I + U) (I - U) (a^i^s) (0+ , y ; T) = 0, (4.8b) 
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where we have used Z2 symmetry Kr{x, —y;T) = K^{—x^y]T). Equations (|4.8a|) and 
()4.8bp are precisely the same boundary conditions for a free particle at the fixed point of 
boundary RG flow. 

To summarize, we have shown that one-particle quantum mechanics on a star graph 
with edge-independent bulk interaction at a fixed point of boundary RG flow is formulated 
into the following configuration space path integral 

l-x{T)=x 

]fC(x, y; T) = V U" / exp(i5[x(t)]), U G ?7(iV), = I, (4.9) 

„=0,1 A(0)=(-l)"y 

where ^[^(t)] is a generic parity invariant one-particle action in the Lagrangian formulation 
and given by 5'[x(t)] = dt L{x{t),x{t)). The most important Lagrangians for physical 
applications will be as follows: 

'U^{t)-V{x{t)), 
L{x{t),x{t)) = I (4.10) 
-x\t) + eA{x{t))±{t), 

where V is a, parity invariant external potential, V{—x) = V{x), and A is an external vec- 

P 

tor potential that transforms under the parity as —A, with e being a gauge coupling 
constant. Examples of such external potentials are harmonic potential V{x) = (l/4)a;^a;^, 
Poschl- Teller potential of cosh type V{x) = —gj cosh^(x), and so on. (Note that in this 
argument bulk theory does not necessarily lie on a fixed point of RG flow.) Scalar kernel 
for these parity invariant regular potentials can be continuous and smooth at the origin 
such that the boundary conditions for the total Feynman kernel K is the same as ()4.8ap 
and ()4.8bp . On the other hand, if one consider singular potentials such as Coulomb po- 
tential V{x) = Ze^/|x|, Poschl- Teller potential of sinh type V{x) = g/ sinh^(x) or inverse 
square potential V{x) = gjx^ (namely, conformal mechanics on star graph), the boundary 
conditions for IC must be modified. Nevertheless, the weight factors are still given by an 
A^-dimensional unitary representation of Z2 as long as scalar kernel itself and its boundary 
conditions at the origin are both invariant under the parity P . 

Time-reversal symmetry. Before closing this Section let us briefly discuss how symme- 
tries constrain the weight factors {W(/), W(P)}. For the sake of simplicity, we will restrict 
ourselves to the symmetry under the time-reversal 7" : T 1— )• — T, which acts on the bulk 
wave function as T : '0(x,T) 1— t- '^\)*'{x^ As is well-known |40j . time-reversal symmetry 

can be broken by spatial boundary even if the bulk theory is invariant under T. We would 
like to classify T-invariant subfamily of weight factors from the viewpoint of representation 
theory of Z2. 

Since the time evolution is given by ipi^^ T) = f^dy K{x, y; T)ip{y), the system becomes 
T-invariant if and only if the following condition holds 

K*{x,y;-T)=K{x,y;T). (4.11) 

Notice that for a unitary kernel the condition ()4.1ip is equivalent to K7'{x, y; T) = IfC(y, x; T). 
Hence, if the bulk scalar kernel is invariant under T, K^{x,y; —T) = K^{x,y;T), the 
condition (|4.1ip to the ansatz ()4.2p leads to the following additional constraints for the 
weight factors 

W*(I) = W(/), W*(P) = W(P), (4.12) 
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which are equivalent to W^(/) = W(/) = W-^{I) and W^(P) = W(P) = W-^{P) for 
hermitian unitary matrices. Thus, for time-reversal invariant systems, the weight factors 
{W(I),W(P)} must be an A^-dimensional orthogonal matrix representation of the cyclic 
group 7j2', that is, it must be a map W : Z2 — )• 0{N). Such matrix representation is given 

by 

W(/) = I, W(P) = O, (4.13) 

where O G 0{N) is a symmetric orthogonal matrix that satisfies = I. Notice that 
()4.13p coincides with the previous results discussed in the case N = 2 |40j. Other sym- 
metry constraints such as invariance under permutation of different edges can be studied 
analogously. 



5 Conclusions and discussions 



In this paper we studied path integral description of one-particle quantum mechanics on 
a star graph and proposed a weight factor formulation of boundary conditions. Our pro- 
posal is based on the folding trick, where in the folded theory wave functions become 
A^-component vector-valued functions. In Section [2] we carefully studied free particle case 
and showed that U{N) family of boundary conditions can be encoded into the momentum 
dependent weight factors {I, in the phase space path integral, where S(j») G U{N) is 

a hermitian analytic unitary S-matrix that satisfies E>{p)S{—p) = I. For the price of weight 
factor formulation, however, we have to include non-standard contour integral with respect 
to the initial momentum p{0) = p in order to pick up the bound state contributions. In 
the subsequent Section we showed that, when the theory lies on a fixed point of boundary 
RG flow, the path integral can be simplified to the configuration space path integral with 
the momentum independent weight factors {1,11}, where U G U{N) is a hermitian unitary 
matrix that satisfies = I. In Section [5] we generalized to the interacting case by consider- 
ing edge-independent parity invariant bulk Hamiltonian and showed that, without referring 
to the boundary conditions, the momentum independent weight factors are generally given 
by an iV-dimensional unitary representation of the cyclic group Z2; that is, scale- invariant 
boundary conditions have emerged only through the Z2 structure of bulk geometry. If 
the system is time-reversal invariant, the weight factors are reduced to an A^-dimensional 
orthogonal matrix representation of Z2. 

We believe that our results presented in this paper give a major step toward the un- 
derstanding of path integral formulation on quantum graph. We also believe that our 
results provide the foundation for worldline formalism [H] on star graphs, which is the 
first-quantization approach to perturbative quantum field theory. Indeed, by using the re- 
sults ()4.10p with ()4.1ip it is straightforward to apply the techniques developed in [8]. For 
example, the ground state energy (Casimir energy) per unit volume at the position {x,j) 
for relativistic T-invariant free massless real scalar field theory on M-^ X (star graph) at 
a fixed point of boundary RG flow is evaluated by the following Wick-rotated {iT — t- T) 
configuration space path integral 

1 f°°dT 1 rx(T)=x f I 



^Casimir(x,j) = -^'^jj (47rT)^/2 / (q) "^^^^^ ^xp ( - / dr^x'^ir] 



1 r(^) ^ 

2(4vrx2)(i5+i)/2^JJ' 



(5.1) 
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where we have subtracted n = term, which is the contribution from the whole hne, and 
the factor 1/{4ttT)^/'^ is the one- loop contribution on M^. Equation ()5.ip coincides with 
the previous result |45j obtained in the operator formalism when D = 1. Similarly, one 
can compute n-point Green's function, one-loop effective action and Schwinger effect for 
relativistic scalar field theory (with edge-independent background gauge field) on a star 
graph. 

We are left with a number of questions, however. One of the main questions is the 
following: Is it possible to derive the weight factors {I, S(p)} as the momentum dependent N- 
dimensional unitary representation of the cyclic group Z2? One may expect that, by allowing 
the weight factors to depend on the initial momentum p{0) = p (that is, weight factors in 
the phase space), the initial condition, composition rule and unitarity for the momentum 
representation of the total Feynman kernel ]K(a:, y;T \ p) will lead to the unique solution to 
the hermitian analytic unitary S-matrix §(p) G U{N) that satisfies S(p)§(— p) = I. However, 
the S-matrix obviously depends on the explicit form of the bulk Hamiltonian H such that in 
order to determine §(p) we need an additional information about bulk interaction. So how 
should we include such information without solving the bulk Schrodinger equation nor using 
the boundary conditions? One of candidates is to impose an additional condition that the 
bulk system is semiclassical and one-loop exact (or WKB exact), which is the case of free 
particle discussed in Section [21 Another important issue is that we have to determine the 
integration contour C from some algebraic relations. In our approach, determination of the 
S-matrix singularity structure is mandatory. Resolutions to these problems are welcome. 
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A Proof of or t honor mality and completeness 



In this Appendix we discuss computational details for the proof of orthonormality and 
completeness relations of the energy eigenfunctions (|2.8p and ()2.10p . Since the relations 
(|2.11bp and (|2.11cp are easily derived, we will concentrate ourselves to the proof of (|2.11ap 
and (j2.12p . We note that essentially the same computations presented in this Appendix can 
also be found in [10113^ . 

Let us first show the orthonormality relation ()2.11ap . Substituting the solution (|2.8p 
into the left hand side of (|2.11ap we get 



dx [V^"(x;p)]^ • V^(x;g) 



sab 



2Tr6{p - q)+ 2tt6{p + q) 



ipL"- + 1 



+ V 
-V 



1 



p — q \ ipL"- 
I I ipL°- + 1 



ipL"- 
ipL"" + 1 iql"" - 1 
1 igL° + 1 
iqL" - 1 



1 



p + q \ipL'^ — 1 iqL"- + 1 



(A.l) 
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where we have used the relation (^")^ • ^ = 6°"^ and the integration formula 

poo ■ • 

/ dxe^^P^ = ± lim = ±V- + tt5{p). (A.2) 

Jo e^o+p±ie P 

Here V stands for the principal value. It is easy to check that the last two lines in (jA.ip 
cancel each other. Noting that 5{p + q) = for g > 0, we obtain the equation (|'2.11a|) . 
Let us next show the completeness relation ()2.12p . A straightforward calculation gives 

X; r^r{x;p) ■ [r{y;p)y = ir ge^^^-^) + T |s(p)e^^'(-+^) 

= 6{x-y)I- -^e-(^+^)/^>^ (A.3) 

0<L"<oo 

where x,y > 0. In the first line we have used = IP"; Z^o^i = I and the definition 

()2.9p . and in the second line the contour integral with the contour along the real p-axis 
enclosing the upper half complex p-plain. Since the second term of ()A.3p is nothing but 
the opposite sign of the contribution from the bound states, X]o<L''<oo V'sC^^) " [V'sl?/)]^; 
we arrive at the completeness relation ()2.12p . Now it is obvious that, as noted in |45] . the 
condition for the absence of bound states is simply expressed as follows: 

r ^ S(p)e^P^ = for X > 0. (A.4) 
J-oo 2vr 
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